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ABSTRACT 

The dynamics of magnetic field decay with Hall drift is investigated. Assuming that 
axisymmetric magnetic fields are located in a spherical crust with uniform conductivity 
and electron number density, long-term evolution is calculated up to Ohmic dissipa- 
tion. The nonlinear coupling between poloidal and toroidal components is explored 
in terms of their energies and helicity. Nonlinear oscillation by the drift in strongly 
magnetized regimes is clear only around the equipartition between two components. 
Significant energy is transferred to the poloidal component when the toroidal compo- 
nent initially dominates. However, the reverse is not true. Once the toroidal field is 
less dominant, it quickly decouples due to a larger damping rate. The polar field at 
the surface is highly distorted from the initial dipole during the Hall drift timescale, 
but returns to the initial dipole in a longer dissipation timescale, since it is the least 
damped one. 
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1 INTRODUCTION 

The recent discovery of soft gamma repeaters (SGRs) with weak fields ( Rea et al. 20ld ) has again raised problems related 
to magnetic field evolution in isolated neutron stars. The activity of magnetars with S GRs and anomalous X-ray pulsar s 



(AXPs) has been believed to be powered by the decay of ultrastrong magnetic fields ( Thompson fe Duncan 19951 . .1996) 
Their surface dipole field, which is observationally inferred from the spin period and its time derivative, typically exceeds 
electron quant um magnetic fie ld Bq = 4.4 x 10 13 G. The field strength of SGR0418+5729 is, however, relatively weak at 
< 7.5 x 10 12 G (|Rea et al.ll20ld ). The critical boundary between magnetars and radio pulsars has thus become less clear, and 
magnetars are not sufficiently characterized by their dipole field strength alone. Their activity may be explained by hidden 
magnetic fields such as poloidal components with higher-order multipoles or internal toroidal components. In either case, the 
fie ld strength should be greater than B = 10 14 -10 15 G, since other energy sources are insufficient (see, for example, the review 
bv lMereghettl (|2008h ). 

The importance of H all drift on the magnetic evolution of neutro n stars was pointed out before evidence of the existence of 
magnetars was available 1 Jones|[l988l ; Goldreich fc Reiseneggerll 19921 ). The effect, which depends on the field strength, becomes 
more important in strong reg imes where B > 10 14 G. Furthermore, Hall drift may in duce instability under certain conditions 
II Rheinhardt fc Geppertj|2002l) . The effect h as been considered in analytic treatments (ICumming et al.ll2004l;lReisenegger et aT 



20071 ). and as plane-parallel slab geometry jVainshtein et al1l2000l ;l Geppert et alj|2003l ; iRheinhardt et al.ll2004l ). These studies 



are useful to understand some aspects of the mechanism, but nonlinear numerical analysis is also required to examine the 
behavior in more realistic stars. Several simulations have also been performed, assumin g that the fields are located in the 
spherical crust region of a neutron s tar jNaito fc Koiimalll994l; Ishalvbkov fc Urpmlll997l ; iHollerbach fc Rudigerll2002i . |2004| ; 
Pons fc Gepperj|2007 ). In particular, Hollerbach fc Riidigei ( 20021 . 20041 ) extensively studied the effect in a crust w i th un iform 
conductivity and density, and subsequently extended the discussion to stratified stellar models. iPons fc Geppertl (|2007l ) also 
calculated the evolution in a realistic stratified model with a thermal history. These numerical calculations are based on 
spectral or quasi-spectral methods, for example expanding angular functions using spherical harmonics. Limitations of such 
an approach are discussed therein. For example, evolution of a purely toroidal field forms a steep gradient that cannot be 
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calculated jPons fe Geppert 2007). A recen t promising approach is using a finite difference scheme to examine the nonlinear 
evolution of the Hall instability in a 2D slab (|Pons fc Gepperti2010] ). Each numerical scheme has advantages and disadvantages, 
so multiple complementary approaches are needed. 

In this paper, we calculate magnetic field evolution using a finite difference scheme to understand nonlinear Hall drift 
dynamics. The model is simplified by assuming an axisymmetric magnetic field located in a crust with uniform conductivity 
and electron number density. This paper is organized as follows. In section 2, the model and its assumptions are described. 
Magnetic field evolution is governed by an induction equation with nonlinear Hall drift term. The relevant boundary conditions 
and initial configurations are also discussed. Sections 3 and 4 give numerical results for two distinct initial configurations. One 
condition is a purely toroidal field in which the poloidal part is always zero if it is exactly zero at the initial state. This is in 
contrast to a purely poloidal initial case, for which the toroidal part is inevitably induced. The evo lu tion of a purely t o roida l 
field is furthermore of interest, due to the similarity to Burgers' equation (e.g., IWhithaml (|l974r )). IVainshtein et al.1 {2000) 
study the problem in a stratified plane-parallel slab, which corresponds to small-scale dynamics much smaller than the stellar 
radius. The results clarify the local mechanism, but our concern is global aspects. How do the results change in a spherical shell 
like the crust? This problem is numerically studied in section 3. Section 4 describes a second study related to the evolution 
of mixed fields with poloidal and toroidal co mponents. Magnetohydrodynamics simulation ( Braithwaite fc Nordlund 20061 : 
Braithwaite fc Spruit 20061 ; Braithwaite 20091 ) shows that dynamically stable configurations are such mixed ones, in which 



poloidal and toroidal field strengths are of the same order. There is little known about the initial configuration of neutron stars 
in particular, the location, topology, and field strength. If both components coexist, their strengths are likely to be similar. 
It is therefore important to explore energy transfer by Hall drift between the components over long-term evolution through 
Ohmic dissipation. Section 5 presents our conclusions. 



2 MODEL AND FORMULATION 



2.1 Equations 

The magnetic field evolution is governed by the induction equation with the Hall term 

„2 



^-B = -V x ( 4— V x B ) + V x 
at \ Aiva 



47ren, 



■B x V x B) 



(1) 



where a is the electric condu ctivity, n e the electron number density, e the charge density, and c the speed of light ( Goldreich fc Reisenegger 
19921 ; Naito fc Koiima 1994 ). There are two typical timescales associated with the first and second terms in eq. {TJ, the Ohmic 
decay timescale AnaL 2 /c 2 and the Hall drift timescale 4-Ken e L 2 /(cBo). Here, L and Bo are typical values for the spatial length 
and magnetic field strength. In general, a and n e depend on the spatial position and the time through the thermal history. In 
this paper, however, for simplicity they are assumed to be constant to explore the dynamics. Thus, the behav ior of this system 
is spe cified by the ratio of two timescales, lZ m = aBo/(ecn e ), which is called the magnetized parameter bv lPons fc Geppert 
(2007). The parameter lZ m is g iven by the initial ma ximum value of the magnetic field. The typical value for a neutron 
star is lZ m =l-10x (_B p/10 13 G) ( Pons fc GeppertJ 12007 ). but this significantly depends on spatial position and temperature 



( Cumming et al.ll2004l ). The overall physical timescale is scaled by Td = 4-7T ar 2 /c 2 , where the stellar radius r s is used for 



normalization. Note that the characteristic decay timescale is Td for node-less field filled in a sphere, but it becomes smaller 
for the field localized in a crust. A magnetic field with lZ m (> 1) is initially given, and the evolution is numerically followed 
until decay at ~ Td- In the case where lZ m 2> 1, which is relevant to strong fields like those of magnetars, the second term on 
the right-hand side of eq. ((TJ) dominates. The advection term is nonlinear, and treating it becomes complicated. 

Magnetic fields with axial symmetry (d/d(f> = 0) are described by two functions, a flux function G describing the poloidal 
magnetic field and a stream function S describing poloidal current flow: 



S ^ 



B = ^(VG x 

where R = r sin 6 is the cylindrical radius in spherical coordinates (r, 

— j = V x B = i(V5 x eV) - ^V(G)e^, 
c K H 

where T){G) is given by 



(2) 



Ampere's equation gives the current density as 



d 1 d 



dr 2 



G. 



or in cylindrical coordinates (R, Z, (j>) by 



(3) 



(4) 



(5) 
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The magnetic field evolution |l} is written in terms of G and S (e.g., Reisenegger et all |2007l )) 



92. = -V{G) + ^-(VG x VS) 

at Td TdR 



(6) 



V x {-L(p{G)VG + SVS)\ 



(7) 



™ = L V{S ) + ^ 

where G, S, and the s patial length are appropriately normalized. As has been pointed out (e.g.. IVainshtein et al. (2000); 
Pons fc Geppertl (|2007l )). the evolutionary equation of a purely toroidal magnetic field is very similar to B urgers' equation, 
which is a simple example of nonlinear propagation with diffusion in one spatial dimension ( Whitham 1974h . Equation (J?} for 
G = is thus reduced to 



dt ~~ T d { dRRdR + 8Z 2 



S + 



2TZ m S dS 
r d R? dZ' 



If the function S depends only on Z, eq. ([8} is exactly Burgers' equation. 

The functions G and S in eq. ([2]) can be expressed by a sum of Legendre polynomials Pi (6): 



G 



!=1 



si(r, t) sin 



g dPi(e) 

89 ' 

} dPi{0) 
89 ' 



The functions gi and si independently evolve in the absence of the Hall term (7Z,-, 
gi and s; with different indices I becomes important with the increase of TZ m . 



decay with the Hall effect h ave been performe d by such an expansion ( Naito fc Koiimal 11994 ; iHollerbach fc Riidiger 



(8) 

(9) 
(10) 



= 0). However, nonlinear coupling among 
Most numer i cal calculations of magnetic 



2002 



2004; |Pons fc Geppertll2007l ). with an exception(Sha lvbkov fc Urpinlll997f ). Angular part of the field is expanded by spherical 
harm onics, but finite difference is used for radial direction in quasi-spectral metho dllNaito fc Koiimal 19941; IPons fc Geppert 
2007). Tchebycheff polynomials are used for radial direction in spectral meth od(|Hollerbach fc Riidigerl 12002 . 2004). The 
limitation of spectral or quasi-spectral methods is also discussed, for example by Pons fc GeppertT i 20071 ). A numerical Gibbs 
oscillation appears when the function evolves to form a steep gradient caused by advection. This is likely to occur in the 
large magnetized parameter lZ m - The finite difference method is used in this paper as an alternative approach. The numerical 
scheme used is a simple stable one, first-order forward time differencing and second-order centered space (FTCS) differencing 
(see, for example, Press et al. 1 19921 )). There are more sophisticated schemes, but it is not easy to apply them to the nonlinear 
Hall drift term, which is most important here. The grid is staggered one with equal spacing. Typical number of grid points 
is 100 x 120. Our numerical results are verified using previous results based on the spectral method in IHollerbach fc Riidiger 
(2002), and the method works well for most parameters. 



2.2 Initial configuration and boundary conditions 

The magnetic field is assumed to be located outside the superconducting core. The crust ranges from ri to the surface r 3 . 
A typical size in neutron star models is (r s — n)/r s ~ 0.1, and depends on the equation of state and the stellar mass. A 
slightly thicker crust model is chosen here, r\/r s = 0.75, because this allows an easier demonstration of the numerical results. 
In actual implementation it is necessary to use more realistic models that include stratified conductivity, number density, and 
so on, but this simple model is useful for understanding the fundamental dynamics. 

The magnetic field cannot penetrate into the core (r < r-y). The condition for the toroidal field is s imply S = at ri. The 



condit i on for the poloidal field m eans that G is a constant at n, which is chosen as G = 0. As discussed in lHollerbach fc Riidiger 



(2002); IPons fc Geppertl (|2007l ) the tangential components of the electric field, E$ and should vanish. They are explicitly 



given by 

The second term on the right-hand side vanishes because B = at n . The first term may be negligible for a large conductivity 
a. Otherwise, some conditions should be imposed to cause the tangential current to vanish. The condition jg — is satisfied if 
the function S rapidly approaches as r goes to r\. The condition j$ = should be imposed on the function G, more precisely 
on T>(G), which is not easily treated. In the numerical simulation, the conditions G = S — are used at ri, assuming a large 
a. The current distributions should be set up at the initial time for the conditions to be satisfied. 

The exterior of the star is assumed to be vacuum. The toroidal field should vanish, so that the boundary condition is 
S = at r„. The vacuum solution of the poloidal field is expressed by a sum of multipole fields, 
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1=1 

where the coefficients a; represent the multipole moments. For example, the dipole moment is a± = fi. The radial component 
B r should be continuous to the exterior vacuum solution across the surface by the condition V • B = 0, while Be may be 
discontinuous if a surface current is allowed. In this paper, neglecting the surface currents, both components are assumed to 
be continuous, and the coefficients ai are calculated by the interior numerical function G at the surface r = r s . The surface 
boundary condition can be expressed by a sum of a; up to lmax =20-30. The results slightly change if the truncation is 
Imax < 10, but do not change even if lmax is further increased. Compared with spectral or quasi-spectral methods, small 
scale structure is not evident in our finite difference method. Initial configuration contains / =1 or 2 component, so that the 
truncation is ju stified. Another approa ch without the multipole expansion is proposed as non-local boundary condition by 
Green's formula! Pons fc Geppertll201ol ). but is not used here. Finally, for the boundary condition on the symmetric axis, = 



and tv are the regularity for the functions G and S, respectively. In other words, G — and S = there. 

The functions G and 5 at the initial time of the numerical simulation should also be subject to the boundary conditions. 
The configuration is easily specified by the functions gi and s; in eqs. © and (|10[) . Initial current distribution is chosen as 
j — at both inner and outer boundaries. One simple solution satisfying this is given by 

a,(r,0) = b ; sin 2 f ^-^A (13) 
^ r,-n J 

where bi is a constant. Both the function s; and its derivative dsi/dr approach as r — > n and r — > r s . For the poloidal field, 
the function gi (r, 0) is numerically given by solving 

d 2 + . f n(r-n) \ .... 

1 gi = - sin , (14) 



dr 2 r 2 J \ r a — n 

with boundary conditions gi(r 1 ,0) = and gi(r 3 ,0) = a;. The source term in eq. p4[) comes from which is chosen to be 
localized near the geometrical center r w (n + r 3 )/2. The function gi(r, 0) is smoothly connected with the multipole solution 
ai/r l at the surface r s . 



2.3 Energy and helicity 

In previous works (e.g., Hollerbach fc Riidigei ( 20021 . 2004 )). coefficients of Legendre polynomials were utilized to study the 
system dynamics, and are useful to understand, for example, energy transfer among different wavelengths at a given time. 
In order to examine the whole dynamics, it is necessary to calculate the coefficients for many times, because they are time- 
dependent. Our concern in this paper, however, is global aspects, so that energy E(t) and helicity H(t) are used to represent 
the magnetic fields. These are indicators of the field strength and twisted structure. Integrating over the entire space, they 
are given by 



E 



8tt 



(B ■ B)dV, 



H = {A- B)dV, 



(15) 



(16) 



where dV is the three-dimensional volume element, and A is a vector potential for B. The energy is divided to poloidal Ep 
and toroidal Et parts, E = Ep + Ep, by which the energy transfer between them is examined. From eq. ([2}, the explicit 
forms are given by functions G and S: 



E P 



E T = 



1 



H = 2 



GS 



dV. 



8tt J \RJ ' J R 2 

From cqs. (O and $J$, the evolution of magnetic energy is given by 



(17) 



(18) 



This is nothing but Oh mic dissipation. The Hall drift does not concern the energy dissipation (jGoldreich fc Reiseneg ger 1992; 
Naito fc Koiimalll994r h but may affect the dissipation rate by modifying the current distribution. Similarly, the evolution of 
magnetic helicity H is 



H = - 2c(E ■ B)dV = 



2c 



(j ■ B)dV. 



(19) 



© 0000 RAS, MNRAS 000, 000-000 



Magnetic field decay with Hall drift 5 



In this way, both energy and helicity decay through Ohmic dissipation. Note that E and its dissipation term, the right-hand 
side of eq. (|18[1 . have a definite sign, but H and its dissipation term, the right-hand side of eq. (|19[) do not. The energy should 
monotonically decrease, but evolution of the helicity is unknown without numerical integration. In the numerical calculations, 
the energy balance equation (|18jl is used to check the accuracy. The relative errors to the initial magnetic energy are less than 
10 -2 . They decrease with increase of the grid number. Relatively large numerical error is produced in the simulation with 
large magnetized parameter lZ m - It is also prominent in the initial evolution, in which steep structure is formed. Functions in 
the late phase are rather smooth, so that the errors do not accumulate so much, even if numerical integration continues to a 
longer time. 



3 EVOLUTION OF PURELY TOROIDAL FIELDS 

This section discusses the evolution of purely toroidal magnetic fields. The system is similar to that of Burgers' equation, as 
discussed in Section 2. The initial configuration of S is specified by a single component si(r, 0) in eq. (|13|) . The three models 
considered here are characterized by I = 1 with b 1 > in Model A, I = 2 with 62 > in Model B, and I = 2 with 62 < 
in Model C. The toroidal magnetic field B^ in Model C has the same amplitude as Model B, but opposite direction. These 
initial configurations are shown in the top to bottom panels of Fig. 1. The coefficient b\ or &2 is chosen for the maximum of 
\B$\ to be equal to 1, since the magnetic fields in the numerical calculation are normalized by the maximum value at the 
initial state. The magnetized parameter for these models is lZ m = 100. 

The top row shows the evolution of £ = 1 with 61 > (Model A). Snapshots of the configuration are shown at the 
initial state (left), at the Hall drift timescale ~ 6 x 10 -4 Td (center) and at the Ohmic d ecay timescale ~ 2 x 10~ 3 rd (right). 
Note that the amplitude of si decays by exp( — 166t/rd) in the absence of the Hall term (jHollerbach fc Riidigei 2002h . so the 



typical decay time is not Td but 6 x 10 _3 rd. The difference mainly comes from the choice of normalization length. The stellar 
radius is used in this paper to compare previous works, but more appropriate one is the crust size. Thus, Td is rather large to 
characterize the actual Ohmic decay. The maximum of S at t = is located at r ~ (ri + r s )/2 and 9 = 7r/2, and moves in the 
negative ^-direction, as shown in close-ups in Fig.l. The drift stops at the outer boundary in several times the drift timescale. 
It is clear that the function S at this time contains higher multipole components s; (r, t) in addition to the initial s\ (r, t) when 
it is expanded as in eq. ()10[) . Subsequently the shape is nearly fixed, but overall amplitude gradually decreases with a longer 
Ohmic decay timescale. Some numerical calculations were performed longer until ~ Td, but the behavior after 6 x W~ 2 Td is 
well described by a simple exponential decay. The results are limited to the early phase since subsequent evolution is easily 
inferred from the extrapolation. 

The middle row shows the evolution of Model B (62 > 0), and the bottom row shows the evolution of Model C (62 < 0). 
The evolution of the magnetic configuration is quite different according to the initial sign of 62- The positive and negative 
regions of S 'collide' at the equatorial plane in Model B, while they 'repulse' each other in Model C. This occurs because the 
drift is negative in the Z-direction for S > 0, but positive in the Z-direction for S < 0. It is also important to note that the 
shape predominantly moves in not the ^-direction, but the Z-direction (see eq. ©)■ The motion in S is clear in close-ups of 
steep structure in Fig.l. In any model, the moved shape decays in a longer decay timescale. 

These results remarkably show the nature of the Hall drift, where a different initial sign for B^ leads to a different fate. 
This leads to an interesting question: Do these different configurations lead to significantly different dissipation rates? If so, 
the direction of B^ would have a significant effect. Figure 2 shows magnetic energy time evolutions for three models. The 
difference between Model B and Model C is relatively minor. A more marked difference is shown for the initial multipole. That 
is, the energies for Model B and C decay slightly faster than that of Model A. The current is swept to the outer boundary 
in Model C, but the outer boundary is replaced by the equatorial plane in Model B. In Fig. 2, a curve with exp(— 332t/rd), 
corresponding to the energy decay in the absence of the Hall drift, is also plotted for comparison. 

Figure 3 shows the dependence of the magnetic parameter lZ m . The initial configuration is the same as that of Model B, 
the model of I = 2 with &2 > 0. The same figure also shows the magnetic energy for free decay. The damping is clearly strong 
for large lZ m . The Ohmic decay is enhanced in the presence of the Hall drift, that is, current is swept into a certain region, 
where it is effectively dissipated. 



4 EVOLUTION OF MIXED FIELDS 

In this section, field evolution is numerically studied for a mixed magnetic configuration consisting of poloidal and toroidal 
fields. The initial configuration is given solely by the I — 1 component for both fields, namely, eq. (|13[) for si and eq. (|14[) for 
g\. The maximum of each field is chosen as the same amplitude, and the magnetized parameter is TZ m = 100. Figure 4 shows 
snapshots of the evolving fields at representative times. The color contour represents the function S of the toroidal field, and 
lines denote the contour of the magnetic flux function G of the poloidal field. 

Oscillatory behavior is clearly evident in G. Initially, the function decreases with the increase in cylindrical distance, and 
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Figure 1. Snapshots of time evolution for the function S, Model A (b\ > 0) at top, Model B (62 > 0) in the middle, and Model 
C (62 < 0) at bottom. Column times are t/rd = 0, 6 X 10 -4 ,2 X 10 -3 from left to right. Whole spherical regions are shown by color 
contour in the left nine panels, while their close-ups of the steep region are shown by contour line with increment of 0.1 in the right 
nine panels. The region is limited to 0.5 < R/r a < 1, -1 < Z/r s < in Model A, 0.5 < R/r s < 1, -0.5 < Z/r s < 0.5 in Model B, and 
0.5 < R/r s < 1, -1 < Z/r s < in Model C. 
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Figure 2. Normalized magnetic energy as a function of time t/rd for three models with magnetization parameter TZ m = 100. Free decay 
curve is also plotted for comparison. 

the maximum is located on the equator 8 = tt/2. The maximum moves 'upward' in the meridian plane, toward 8 < tt/2, 
until t/rd ~ 1-4 x 10 -3 (second panel). It then changes direction and goes 'downward,' passing through the equator at 
t/rd ~ 3.2 x 10~ 3 (third panel) and reaching a minimum at t/rd ~ 5.2 x 10 -3 (fourth panel), before returning to the initial 
position at t/rd ~ 7.8 x 10 -3 (fifth panel). During this cycle, the field strength decreases. 

The function S is also oscillatory. The initial configuration contains only the I = 1 component in the angular part 
(5* oc sin 2 #), which is symmetric with respect to 8 — tt/2. The state at t/rd ~ 1.4 x 10 -3 (second panel) markedly differs 
from the initial state. The configuration is no longer symmetric, and higher multipoles can be seen. The field strength itself 
is weak around this time. At t/rd ~ 3.2 x 10~ 3 (third panel), the configuration again becomes symmetric like the initial 
state, but the sign of S is reversed. The I = 1 component is dominated there. After the direction of B^(= S/(r sin 8)) again 
changes, the configuration returns to the initial one at t/rd ~ 7.8 x 1CP 3 (fifth panel). The directional change occurs around 
t/rd ~ 1-4 x 10~ 3 (second panel) and 5.2 x 1CP 3 (fourth panel), which correspond to a local minimum of toroidal field strength. 
The overall toroidal field strength also decreases during this cycle. 

Figure 5 clearly shows the oscillatory behavior of the magnetic energy, which is divided into poloidal and toroidal parts, 
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Figure 3. Normalized magnetic energy as a function of time t/rd for TZm = 20, 100, 500. Free decay curve is also plotted for comparison. 




Figure 4. Snapshots of time evolution for functions G and S. Contour lines outwardly represent the level of G for 0.02 X n X (Bqt^ ), n = 
1,2, ■ • ■. Color contour represents S normalized by Bor s . Note that different color scales are used, since S becomes very small at the 
bounces in the second and fourth panels. Those panels use the color scale on the left; the others, the scale on the right. 



Ep and Et, respectively. The magnitude of the amplitudes at t = are nearly the same. □ The curve of toroidal energy Et 
represents a damped oscillation. Local minima can be seen at t/rj w 1.4 x 10 -3 and 5.2 x 10 -3 , corresponding to the time of 
the second and fourth panels in Fig. 4. The configurations of the third and fifth panels in Fig. 4 are those of the local maxima 
at t/rd ~ 3.2 x 10~ 3 and 7.8 x 10~ 3 . There is remarkable energy transfer between the toroidal and poloidal parts during the 
initial phase. Initially Ep increases until t/ra « 1.4 x 10 -3 , although the total energy E aum = Ep + Et decreases. The sum 
always decreases due to Ohmic decay (see eq. HE}). The initial rapid decay of Et is thus partially due to this transfer. Energy 
is subsequently transferred between the two components in turn, but the behavior becomes less clear. The magnetic field 
decays on the timescale t/rd ~ 1CP 2 , so the coupling becomes weak. The oscillation period gradually becomes longer, since 
the drift timescale increases. Figure 5 also shows the evolution of magnetic helicity H, which exhibits oscillatory damping 
with Ohmic decay timescale t/rd » 10 -2 . The change of sign in H denotes an inversion of the toroidal field B^, which occurs 
around the local minima of Et- 

The coefficients a; of the multi-moments in eq. (|12[) describe the exterior poloidal field. The right panel of Fig. 5 shows 
the evolution of a few of the lowest values. Higher multipoles are induced until t/rd ~ 1.4 x 10 -3 , the configuration shown in 



Values are not exactly the same, because the maximum amplitude of each field is fixed at same value, the distributions are slightly 
different. 
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Figure 5. Time evolution of energy E, helicity H (left panel), and multipole moments a; at the surface (right panel). Energy is 
normalized by the initial total energy, and helicity by the initial value. Coefficient a; of the l-th moment is normalized by the initial 
dipole value air l s ~ . 



the second panel of Fig. 4. The poloidal field is no longer dipole at this time. Interestingly, the coefficient ai is significantly 
decreased, as compared with the initial value, although a large amount of energy is stored in the poloidal part, as shown in 
the left panel of Fig. 5. Dipole field strength at this time is not a good indicator of overall magnetic energy. The poloidal 
field returns back to the dipole at t/rd ~ 3.2 x 10 (the time of the third panel in Fig. 4), at which higher multipoles are 
temporarily zero. After sinusoidal oscillation, higher multipoles decay rather rapidly, so the polar dipole remains for the later 
time t/rd > 10~ 2 . The behavior after 10~ 2 Td is rather simple. Magnetic energy is significantly dissipated, so that the Hall 
drift becomes less important. Subsequent evolution in the late phase is described by a free decay due to Ohmic dissipation, 
and is omitted here. 

Figures 6-9 compare evolutions of energy, helicity, and polar multipole coefficients for four models. The initial magnetic 
configuration is the same as that of Figs. 4 and 5, but the strength and magnetized parameter 1Z m vary. The ratio of poloidal 
to toroidal fields is 1/4, or 4 in amplitude, so 16(= 4 2 ) or 1/16 in energy, and lZ m is 7Z m = 20, or 100. The right panel shows 
energy plotted on a logarithmic scale, log 10 (_E), and the left helicity H and a few lowest multipole coefficients ai. Figure 6 
shows the results of a model with 1Z m = 20, in which the energy is initially dominated by the poloidal part. Polar multipole 
components are induced through the coupling to the toroidal part, but the dipole field is barely affected. The decay curve 
is well described by exp(— 55t/rd), which coincides with free decay of the dipole. The energy is always dominated by the 
polar dipole, and the evolution is described by exp( — llOt/rd), a line in the logarithmic scale. The curve of toroidal energy 
is nonlinear on the logarithmic scale, but is oscillatory. The minor component is highly affected by the poloidal one, but is 
almost neglected in the dynamics of the whole system. No consequence of the toroidal field comes from the fast decay rate 
exp(— IQQt/rd) in the amplitude as seen in Section 4, or from the initial small strength. 

Figure 7 shows the results of an evolution in which the toroidal component is dominated at the initial state. Compared 
with the decay curve of total energy in Fig. 6, the damping is much faster. Most of the magnetic energy is initially stored 
in the toroidal part, but rapidly decays. The poloidal component decays rather slowly, and dominates at the later time. The 
higher multipole moments are induced and have larger amplitudes than those of Fig. 6. The overall dipole decay curve is 
nonetheless very similar to that of the large poloidal case. 

Figure 8 shows the results for large a magnetization parameter 7Z m = 100 and an initially large poloidal field. The 
evolution of the total energy is very similar to that of Fig. 6, although the oscillatory behavior in the toroidal energy is 
evident. The timescale determined by the Hall drift is approximately 1/5 that of Fig. 6, in which TZm = 20. The oscillatory 
behavior is also clear in the coefficients a; and helicity H. However, except for the initial wavy structure seen in the Hall drift 
timescale, the decay of the dipole field is similar to that of Fig. 6. The dominant component, the polar dipole, is thus not 
affected by the toroidal field, which never plays an important role because of its rapid decay. 

Figure 9 shows the results for a large magnetization parameter lZ m = 100 and an initially large toroidal field. In this case, 



© 0000 RAS, MNRAS 000, 000-000 



Magnetic field decay with Hall drift 9 




Figure 6. Time evolution of energy in logarithmic scale log 10 (E) (left panel) and helicity H and multipole moments a; at the surface 
(right panel) for an initially large poloidal field with TZm = 20. These normalizations are the same as in Fig. 5. 




Figure 7. Same as Fig. 6, but for an initially large toroidal field with TZm = 20. 



a large amount of energy is transferred from the toroidal to poloidal components. Higher multipoles in the poloidal field are 
induced around t/rj ~ 10~ 3 . The configuration of the poloidal field at this time is significantly distorted as in the second panel 
in Fig. 4. The toroidal energy becomes much smaller than the poloidal energy around t/ra ~ 10~ 2 , and both components are 
subsequently decoupled. After that, the dipole field evolution is determined by free decay. Despite the disorder at t/ra « 10 -3 , 
the amplitude of the free decay phase, for example at t/ra = 2 x 10 -2 , does not differ so much from that of Figs. 6-8. 



5 SUMMARY AND DISCUSSION 

Numerical simulations demonstrate how the Hall drift changes the current and magnetic configuration. In a purely toroidal 
evolution, Ohmic dissipation is enhanced by accumulated currents elsewhere. The spatial location depends on the initial data, 
but the energy dissipation rate does not so significantly depend on the accumulation position in our uniform conductivity 
model. In a realistic case the conductivity will decrease with the radius, so most of the magnetic energy may be effectively 
dissipated near the surface. The results become more sensitive to the initial condition. Energy is, in principle, transferred 
between poloidal and toroidal components if both are initially involved. It is important for understanding the evolution that 
the free decay rate of the polar dipole is the least. When the polar dipole dominates, it is rarely affected by the toroidal 
field, which rapidly decays. On the other hand, a significant amount of energy is transferred to the poloidal field until almost 
equipartition when the toroidal component is initially dominant. Global nonlinear coupling is manifest in the Hall drift 
timescale only when the corresponding energies are of the same order. Moreover, the poloidal field at the surface and the 
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Figure 8. Same as Fig. 6, but for an initially large poloidal field with lZ m = 100. 




Figure 9. Same as Fig. 6, but for an initially large toroidal field with TZ m = 100. 



exterior is highly distorted, and is no longer described by a pure dipole field. During this phase, the dipole is not a good 
indicator of overall field strength. In a longer Ohmic dissipation timescale the toroidal field decays rapidly, so the coupling 
vanishes. The polar dipole eventually survives. 

The early evo lution highly depends on the choice of the initial data, i.e, the configuration and strength. As discussed by 



Reisenegger et al.l (|2007r ). stationary conditions in the presence of Hall drift for an axisymmetric magnetic field means that 



the isosurface of S coincides with that of G; in other words, S = S(G), for which the coupling term VG x VS 1 vanishes. 
This condition is not satisfied in the whole shell region. The outer boundary condition requires that the toroidal component 
is concentrated in the interior (5 = 0), whereas the poloidal one may leak out to the exterior, meaning that G / at the 
surface. The interior poloidal field is described by the value at the surface if the topology is simply connected. The condition 
S(G) therefore means that S = everywhere, irrespective of G. Our numerical calculation starts with an initial configuration 
different from the 'Hall equilibrium', and shows the behavior to the state S — 0. It may be necessary as a next step to study 
plausible initial configurations, since little known about them. 

The timescale is less accurate in our simplified model with uniform density and conductivity. Realistic models are neces- 
sary, but the distributions depend on equation of state and neutron star mass. Assuming that our uniform model is obtained 
as a result of spatial average of a certain model with stratified number density and conductivity, the timescale is estimated. 
The overall normalization constant is given by = 47tot 2 /c 2 « 4.4 x 10 9 (<t/10 25 s -1 ) (r s /10km) 2 years, where a is an 
averaged value of the conductivity. The decay timescale of the dipole is 2 x lQ~ 2 Td ~ 9 x 10 7 years, and that of the toroidal one 
is 6 x 10 _3 Td ~ 3 x 10 7 years. These numbers are slightly larger, since the crust size in our present model is thick, L = r s /4. 
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The timescale is proportional to L 2 , so that the actual values may be smaller by a factor of ~ 10 _2 -10 _1 . A typical decay 

timescale is around 10 5 -10 7 years. 

The characteristic age of the low-field magnetar SGR0418+5729 is more than 2.4 x 10 years jRea et alJlioich . The 
dipole field (< 7.5 x 10 12 G) may decay or survive within this period. However, the internal toroidal field is likely to dissipate 
more quickly. It is therefore difficult to understand why the toroidal field has a much larger field strength, > 10 13 -10 14 G, 
which is required for the activity. One possible explanation is that current age estimates are inaccurate. The characteristic 
age is normally estimated by assuming a constant dipole magnetic field. The Hall drift significantly affects the surface value, 
especially at the initial epoch with strong field. As demonstrated in the second panel of Fig. 4, the surface field is highly 
distorted from the pure dipole. SGR0418+5729 may correspond to a young phase of oscillatory evolution in which the surface 
dipole temporarily decreases, but there is a strong internal toroidal component. 
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